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Abstract 



^ ' We calculate the critical coupling 4/^^ and critical exponent /3 for 

Q . the order parameter in SU(2) lattice gauge theory by applying of the 

finite size scaling technique and the method proposed by Kouvel and 
Fisher for analysis of experimental data. In contrast to the standard 
finite size scaling approach, this method allows to determine simul- 
taneously both (7c /4 and /? as two parameters of the linear fit to the 
Monte-Carlo data. 
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O I Lattice calculations by Monte-Carlo (MC) method have so far been the 

only method for studying the temperature phase transition in gauge theories 

■ from the first principles. By their nature, MC calculations are carried out on 
the finite lattices. This is a most serious drawback of a computer simulations, 

. because no finite system can exhibit a true phase transition. Nevertheless, 

finite systems remind of phase transitions and systematic studies of these 
Qj \ pseudo-transitions as function of system size may reveal information about 

■ the phase transition in the thermodynamic limit. One way to do this partic- 
ularly in order to evaluate the critical exponents at the theory is to use finite 

^ . size scaling (FSS). 

FSS was proposed in statistical physics [1, 2] and applied to high precision 
MC data of Ising [3, 4] and other models. The validity of this method for 
gauge theories has been demonstrated through the investigation of critical 
properties of SU(2) lattice gauge theories at finite temperature [5-9]. 

In this letter we propose the method of analysis of MC data which is 
the combination of FSS technique and so-called Kouvel-Fisher method for 
analysis of experimental data in solid state physics [11]. We consider MC 
data for SU(2) order parameter on the lattices 8^ x 4, 12^ x 4 and 18^ x 4. 
In contrast to the standard FSS approach, our method allows to determine 
the critical coupling i/g"^ and critical exponent P as two parameters of the 
linear fit. 
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The partition function for SU(2) gauge theory on x N^. lattice is 

n dU,,.e-'^^\ (1) 

where ?7^,,y are the SU(2) matrices on the hnk (//, i/), 5* is the standard Wilson 
action 

SiU) = ^ E (l - IttUp) . (2) 

Up is a product of link matrices around a plaquette. The temperature T (vol- 
ume V) is defined by size of the lattice in the time-like (space-like) direction 
(a is a lattice spacing) 

T = 1/ (iV,a) , V = {N^af . (3) 

The order parameter of the deconfinement phase transition on an infinite 
volume lattice is the expectation value of lattice average of the Polyakov 
loop 

where 

Lr. = ^Tr n Ur,n;0 (5) 
^ r=l 

and f/r,n;0 is the SU(2) link matrix at the point (n, r) in the timc-likc direc- 
tion. On the finite lattice one considers as the order parameter {\L\), since 
in this case (L) is always zero due to system flips between the two ordered 
states. 

According to FSS theory the scaling function for the order parameter on 
the lattice x N^, fix is given by [2] 

m{\L\) = QL{x,g,N^y^), (6) 

where /3, u are the critical exponents for order parameter and correlation 
length correspondingly, gi is the irrelevant scaling field, ~ 1. The scaling 
variable x is defined as 
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Our definition of reduce temperature t differs from one used in [7-10] 

and seems to be more natural because the coupling constant plays just 
the role of a "temperature" of the effective theory which can be obtained 
from the partition function (1) by integration on all link variables except 
those for Polyakov loops. It is clear that this difference disappears in the 
limit 5f^/4 ^ g'^/A . 

At large N^- the scaling function Ql must behave as 

Ql ~ (xf . (9) 

This behaviour may be observed only if x becomes " sufficiently large" . De- 
viations from (9) are in close vicinity of the critical point due to finite size 
rounding and further away from a; = due to correction to scaling. Expand- 
ing Ql in (6) around x — at large N^^ results in 

Nh\L\) = a + ^+[c+^y, (10) 

where a, b, c, d are unknown constants. Dropping the term 0{N^y^) one can 
obtain 

Ar|(|L|) = a + cr, (11) 

what is exactly the approach of Refs. [7-10]. On the other hand at x "suffi- 
ciently large" the x-dependence is drastically changed. In Ref. [6] the scaling 
function has been calculated in the form 



Ql = Ax^ 



[12) 



by the data collapsing method. It has been shown that the correction to 
scaling term in (12) is negligibly small as compared to the leading term on 
the interval 0, 5 < a; < 2. In both cases the SU(2) Monte-Carlo data for the 
order parameter have been used which were computed on x A^^ lattices 
with A^g-=8, 12, 18,26 and — A and have already reported on in Refs. [5- 
10]. The numerical value of the critical exponent /3 has been found which is 
compatible with the corresponding value (3 — 0, 327 of the 3-d Ising model 
[4]. 
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These two procedures have two disadvantages: i) 3 — d Ising model value 
of z/ = 0, 631 has to be used; ii) the very accurate calculation of critical 
coupling is necessary. The latter is determining by the fourth cumulant 
of (|L|) as in Refs. [5, 7] or by the x^-method as proposed in Ref. [8]. 
Even small shift leads to the additional uncertainty in the determination of 
critical exponent (3, because it depends very sensitively on the value inserted 
for ^2/4. 

We propose below the method for the determination of the critical ex- 
ponent /? which does not require to know oi g^/A from the very beginning. 
No information about the critical exponent u is also required. Following the 
Kouvel-Fisher idea which was used for analysis of the experimental data in 
solid state physics [11] we consider the function 

According to the results of Ref. [6] in the region 0, 5 < x < 2 it is easy to 
find from (12) 

+ (14) 



Dropping the correction to scaling term we have in this region with the same 
accuracy as in (11) 



X 



K{x) = -. (15) 



Taking into account the identity 



aq, _,l^-'^aiHL (16) 



dx 4 9(c/V4) 

we have immediately 



K(x) ^-^-JIL m (17) 



Finally, we obtain the simple but informative expression 

(1^1) _ 1 (g' 9l\ 

d{\L\)/d{g^/A) p\a a)- 



(18) 



Eq. (18) means that for the lattices with A^^- fix and various N^j there are 
the intervals of g^/A where MC data for different N„ must arrange around 
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the same straight hne. The inverse slop of this hne represents the critical 
exponent /? and the intersection with 5'^/4 axis yields the critical coupling. 

In order to verify this approach we have used the high precision MC data 
produced by Bielefeld group on the lattices 8^ x 4, 12^ x 4 and 18^ x 4 [5, 
10]. The calculation of derivatives with respect to /A has been made by 
using the parabohc fit for (|L|) (every parabolae includes five MC points). 
This fit does not smooth the statistical fluctuations, but connects the data 
continuously and with continuous first and second derivatives. The least 
minimal x^-fit of Eq. (18) {x^/N = 0, 62) includes 7 MC points in the region 
0.396 < ^74 < 0.426 (2.345 < A/g^ < 2.527) for = 8, 24 MC points in 
the region 0.414 < g^/i < 0.430 (2.325 < 4/^^ ^ 2.415) for = 12 and 16 
MC points in the region 0.427 < g'^/A < 0.433 (2.310 < i/g^ < 2.340) for 

= 18. The result is 

f3 = 0.326 ± 0.003; A/g^ = 2.2988 ± 0.0007. (19) 

We see that the critical exponent j3 is in excellent agreement with 3 — d Ising 
model value and the inverse critical coupling A/g"^ coincides with that value 
which has been calculated in Refs. [7-10] by other method. 

In conclusion the following comments are in order. We have shown that 
the Kouvel-Fisher modification of FSS gives the extremely useful method for 
the calculation of critical exponent (3 in SU(2) lattice gauge theory. This 
method allows to investigate the finite lattice system at small but finite 
t — (g^/A—g'^/A) / {gl/^) and large A^^ with infinite volume formulae (18). The 
critical coupling g"^/^ and critical exponent (3 are simultaneously determined 
as two parameters of the linear fit. It is clear that this approach may be 
applied to the calculation of all critical exponents in SU(2) gauge theory and 
to other models which undergo the second-order phase transition. 
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